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Abstract 

We present an exposure-time calculator (ETC) for the Immersion Grating 
Infrared Spectrograph (IGRINS). The signal and noise values are calculated 
by taking into account the telluric background emission and absorption, the 
emission and transmission of the telescope and instrument optics, and the 
dark current and read noise of the infrared detector arrays. For the atmo¬ 
spheric transmission, we apply models based on the amount of precipitable 
water vapor along the line of sight to the target. The ETC produces the 
expected signal-to-noise ratio (S/N) for each resolution element, given the 
exposure-time and number of exposures. In this paper, we compare the sim¬ 
ulated continuum S/N for the early-type star HD 124683 and the late-type 
star GSS 32, and the simulated emission line S/N for the H 2 rovibrational 
transitions from the Iris Nebula NGC 7023 with the observed IGRINS spec¬ 
tra. The simulated S/N from the ETC is overestimated by 40—50% for the 
sample continuum targets. 
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1. INTRODUCTION 

The Immersion Grating Infrared Spectrograph (IGRINS) is a cross dis¬ 
persed near-infrared spectrograph designed using silicon immersion echellc 
grating (e.g., Marsh, Mar & Jaffe, 2007; Wang, Gully-Santiago & Deen et 
al., 2010; Gully-Santiago, Wang & Deen et ah, 2012). This instrument can 
cover the whole H and K bands (1.4 — 2.5 /irn) in single exposure with a 
resolving power of R = 40,000. Two volume phase holographic (VPH) grat¬ 
ings are used as cross-dispersing elements in the H and K bands. IGRINS 
employed two 2048 x 2048 pixel Tcledyne Scientific & Imaging HAWAII- 
2RG detectors with SIDECAR ASIC cryogenic controllers. Figure 1 shows 
the optical design sketch of IGRINS. This instrument has been employed on 
the 2.7 m Harlan J. Smith Telescope at the McDonald Observatory since 
2014 March. Detailed descriptions of IGRINS can be found in Yuk, Jaffe h 
Barnes et ah (2010) and Park, Jaffe & Yuk et ah (2014). 

Exposure-time calculators (ETCs) play an important role in preparing for 
observations by accurately simulating the expected exposure-time to achieve 
a required signal-to-noise ratio (S/N). Examples of ETCs include those for 
MOSFIRE 1 , CRIRES 2 , SOFI 2 , KMOS 2 , HSpot 3 , and the Wide Field Plane¬ 
tary Camera 2 (WFPC2) of the Hubble Space Telescope (e.g., Bernstein et 
ah, 2002; Li Causi, De Marchi & Paresce, 2002). 

We developed the IGRINS ETC by taking into account the telluric back¬ 
ground emission and absorption, the thermal emission and throughput of the 
telescope and instrument optics, and the dark current and read noise of the 
infrared detector arrays. This tool is a part of the IGRINS software which 
includes the IGRINS data reduction pipeline package (Sim, Le & Pak et ah, 
2014). The software was written in the interactive data language Python, 
and the user controls specific calculation modes and input parameters via 
a graphical user interface. In this paper, we present detailed calculation 
algorithms and the results of the IGRINS ETC. 


1 http://www2. keck, hawaii.edu/inst/mosfire/etc.html. 

2 See the ESO exposure-time calculators at http://www.eso.org/observing/etc. 

3 http: //herschel.asdc.asi.it/index.php?page=proposal.html. 
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2. SENSITIVITY SIMULATION METHODS 

2.1. Throughput 

We define the total throughput, T poin t for point sources as 


A°in 1 Atmo Arrirror A r Aph ^optics Alit7q e (t J 

Axtend Atmo AnirroUgr Aph ApticsAje (2) 

where r terms are the throughput values of the average atmosphere (r atmo ), 
telescope mirrors (r m i rror ), immersion grating (r gr ), VPH gratings (r vp h), other 
IGRINS optics (r op tics), slit (r s ii t ), and detector (r qe ). Table 1 shows the es¬ 
timated values of these parameters in equations (1) and (2). We estimate 
that the average throughput value of the atmosphere is 0.77 and 0.83, as¬ 
suming that the precipitable water vapor (PWV) is 4 mm, at the wavelength 
ranges of 1.4—1.9 p m and 1.9—2.5 pm for the H and K bands, respectively. 
The model of the atmosphere transmission is shown in Appendix A. Since 
the point spread function (PSF) on the slit closely resembles a Gaussian, 
we can easily estimate the slit loss for a given image size. We define a slit 
throughput, r s iit, as 
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( 3 ) 


where E 0 is the unobscured flux from the target, and E s i it is the flux passing 
through the slit (Lee & Pak, 2006). If the PSF is a circularly symmetric 
Gaussian function, r s i it can be expressed as the error function: 
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where Y is a half of the slit width and a is a Gaussian width, i.e., a 
= A?/fwhm / (2\/2 In 2), where Aj/fwhm is the measured full width at half 
maximum along the slit length. 


2.2. Background Emissions 

We calculate the background emissions from the Zodiacal light, the Moon 
light, OH emission lines, and thermal emission from the atmosphere. Table 
2 shows the parameters used in the calculations. The number of electrons 
s^ 1 pixel -1 by the background emissions is given by 
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£Zod,Moon,OH — -^'T'extend —T-— ($Zod + ^Moon + 0.1 $oh) (5) 

AAbW-K 

where r ex tend is the total throughput for an extended source, R is the 
spectral resolution (R = 40,000) for IGRINS, A 0 is the central wavelength of 
the H or K filter, AAbw is the band width of the filter, t&zod, d’Moon, and 
4>oh are photon fluxes from the Zodiacal scattered light, the Moon light, and 
OH airglow in units of photon s -1 m -2 arcsec -2 , AVt is the aperture of the 
telescope multiplied by the solid angle on the sky viewed by one pixel (HO 
= 1.00 x 10 -11 m 2 sr pixel -1 ). 

In Table 2, we list the average values of the Zodiacal light 4 at ecliptic lat¬ 
itudes of 0°, 30°, 60°, and 90°. In the simulation, we assume the heliocentric 
ecliptic longitudes at ±90° since the Zodiacal light is not variable at different 
ecliptic longitudes. In the ETC, the average values of the Zodiacal light is 
‘hzod = 55, 41 photon s -1 m~ 2 arcsec -2 for H and K bands, respectively. 

The Moon light is another natural source that is the major contributor 
to the sky background (e.g., Noll, Kausch & Barden et ah, 2012; Jones, Noll 
& Kausch et ah, 2013). We calculate the average values of the Moon light 4 
at the Moon phases of the half Moon and the full Moon for the cases of the 
Moon and target separations at 0°, 30°, 60°, 90°, and 120° (See Figure 2). 
In the ETC, we assume that the target separation is 60° for the cases of the 
half Moon and the full Moon phases. We may consider to apply for all the 
cases of the Moon light in the next versions of the ETC. 

We calculate the $oh value based on OH emission lines which are shown 
in Figure A.11 in Appendix A. We sum up the OH line intensity values of 
whole band. The OH value, used in equation (5) is the summation of OH 
line intensity in H and K bands, 4 >oh = 3390 + 1096 = 4486 photon s -1 
m -2 arcsec -2 . We also consider that these is ratio of photons from OH lines 
that are scattered inside the IGRINS optical box and distributed as scattered 
light across the detector. By using Flat image which is subtracted by Dark 
image (See Figure 3), we calculate the average intensity of pixel values in the 
order regions, and the inter-order regions, then estimate the scattering ratio 
of photon from OH lines inside the optics (See Table 3). The calculated ratio 
value is 0.1 for both H and K bands. 

The background thermal flux in units of electrons s -1 pixel -1 is calculated 


4 See the Sky Model Calculator at http://www.eso.org/observing/etc. 
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by 
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where h is the Planck constant, and B vT is the Planck function of fre¬ 
quency v and temperature T in units of W m -2 Hz -1 sr” 1 . We define the 
total emissivity e as 


£ = 


(1 Patmo) T [ ^-mirror “1“ (1 ^window) /^window] /^ir 


where £ mirror is 0.25, r u 
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2.3. Signal-to-Noise for Continuum 

We observe point sources using the “Nod-on-Slit” mode, where we take 
spectra at two positions along the slit and subtract them from each other 
(Sim, Le & Pak et ah, 2014). This subtracts various sources of background 
while simultaneously gathering signal from the science target during all ex¬ 
posures. 

The continuum signal is 


n e 


Texnhi-Tr, 
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where S'cont is in units of electrons per resolution element, F ZM is the 
zero magnitude flux density in units of W m -2 Hz” 1 , m is the apparent 
magnitude of the continuum source, n exp is the number of exposures, t exp is 
the exposure-time, and r po i nt is the total throughput for the point source. 

The electrons collected by the detector exhibit as a Poisson distribution. 
The noise is calculated from the formula as 


iV C o„t = [Q + ^Cont] 1/2 (9) 

where IVcont is in units of electrons per resolution element, and 

Q ^PsliWexp [^exp (^Zod,Moon,OH T ^dark T ^thermal) T ^read] (10) 

In this formula, p slit is 3.66 pixels in the spatial and spectral resolutions, 
<^dark and d r ead are the dark current and read noise of the infrared detector 
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array, £zod,Moon,OH and ^thermal are given by equations in (5) and (6). For 
more realistic simulations, we should have considered the number of pixels 
in spatial and spectral resolution elements instead of p 2 lit pixels. Since we 
always subtract a background frame, the noise term has a factor of 2 in 
equation (10). 

2-4- Signal and Noise for Emission Lines 

The ETC can simulate the IGRINS observations of emission lines given a 
rest wavelength, flux, doppler shift and width. We use the following defini¬ 
tions for the signal and noise of the emission line (Pak, Jaffe & Stacey et ah, 
2004): 


Shine — AApi X el ^ (,/l.,i — /c) (H) 

i 

N Une = AA P i xe icr(/c,i) y/ n L (1 + l/n c ) (12) 

where /^i is the line flux data samples in units of W m -2 //m _1 , fc is the 
average continuum adjacent to the emission line, AA p i xe i is the wavelength 
range corresponding to the pixel width at line center, <v(/c,i) is standard 
deviation of the continuum, and n l and nc are the sample numbers in the 
emission line and continuum bands, respectively. 

3. SIMULATION RESULTS 

3.1. PHOTOMETRIC MODE 

In the photometric mode, the ETC returns the average S/N for the H 
and K bands, assuming that the target is a point source. In this mode, the 
S/N values are computed for whole spectral H and K bands. The users can 
input the precipitable water vapor (PWV), the exposure-time, the number 
of exposures, the expected K band magnitude, the effective temperature of 
the source, the seeing size, and the Moon light phases. 

3.2. SPECTRAL MODE 

The spectral mode allows the users to view plots of the calculated S/N 
as a function of magnitude or wavelength. Figure 4 shows the plot of S/N 
vs. limiting magnitude for continuum point sources in the H and K bands. 
Figure 5 shows the plots of S/N vs. wavelength in the H and K bands for a 
blackbody source with a temperature of 6000 K, H-magnitude of 12.1, and 
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K-magnitude of 12.0. In this spectral mode, the S/N values are variable 
with wavelength. In this mode, the throughput values of the atmosphere 
and the background signal and noise from OH emission lines are depended 
on wavelength. The model atmosphere is shown in Appendix A. Figures A.9, 
A.10, and A.11 show the atmospheric transmission and OH emission lines in 
the H and K bands. 

4. TEST OBSERVATIONS 

IGRINS had three commissioning runs, in 2014 March, May and July, on 
the 2.7 m Harlan J. Smith Telescope at McDonald Observatory. In order 
to compare our ETC simulation results with real observational results, we 
chose an early-type AOV star, HD 124683, and a late-type K5 star, GSS 32, 
for point source continuum targets, and a reflection nebulae, NGC 7023 for 
an extended line emission target. The data of HD 124683 and GSS 32 were 
taken by using the “Nod-on-Slit” mode on 2014 May 27. The total exposure- 
times are 480 s and 960 s, respectively. NGC 7023 was observed by using the 
“Nod-off-Slit” mode on 2014 July 13. The total exposure-time is 1200 s. 

NGC 7023 is a typical example of a photodissociation region (Lemaire, 
Field & Gerin et al., 1996; Lemaire, Field & Maillard et ah, 1999). The spec¬ 
tra of this target show many strong narrow emission lines from the molecular 
rovibrational transitions of H 2 (e.g., Martini & Sellgren, 1997; Martini, Sell- 
gren & Depoy, 1999). The observed H 2 emission lines from this target are 
good examples to test the simulated emission lines from the ETC. We use 
the AOV standard star, HD 155379, to do the absolute flux calibration for 
NGC 7023 (see Table 4). The flux calibration processes are based on the 
method of Lee & Pak (2006). 

4-1. Signal-to-Noise for Continuum Source 

We compare the S/N of an early-type AOV star, HD 124683, and a late- 
type K5 star, GSS 32, to the S/N calculated using the ETC. We applied the 
same exposure-time, the number of exposures, K-magnitude, and the effective 
temperature of the sources in the ETC. We also applied the measured FWHM 
of the point spread function in the simulations. PWV is assumed to be 2 mm. 
Figure 6 shows S/N of the real observational data and the simulated S/N of 
HD 124683. We estimate that the average S/N of the real observational data 
is 400, and the average value of the simulated S/N is 600. Then, the ratio of 
the simulated S/N and the real observational data S/N is 1.5, corresponding 
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to 50%. Figure 7 shows the S/N values vs. wavelength of GSS 32. The 
average value of the simulated S/N and the real observational data S/N is 
350 and 250, respectively. Then, the ratio is 1.4. The simulated S/N of 
GSS 32 is overestimated by 40%. 

The observed spectra of the stars are curved, while the simulated spectra 
are flat because we did not apply the blaze function to display the grating 
efficiency of echelle orders in the simulations. In the simulations, we focus 
on the calculations of the expected average S/N values, and simply compare 
that to the S/N values of the real observation data. The actual numbers 
would be different by a factor of 1.5 depending on the blaze angles. 

4-2. Signal-to-Noise for Emission Line 

We compared the simulated spectrum of NGC 7023 with those from the 
ETC. We applied the same exposure-time, the number of exposures, K- 
magnitude, and the effective temperature in the ETC. We chose a sample of 
the brightest emission line H 2 1-0 S(l) 2.12183 pm and the fainter emission 
line H 2 7-5 0(5) 2.02200 pm. In the ETC, we used the intensity values of the 
emission lines from Martini & Sellgren (1997). The flux values which applied 
to the simulations of the IGRINS ETC are 12.8 x 10~ 18 and 0.81 x 10“ 18 
W rn~ 2 for H 2 1-0 S(l) and H 2 7-5 0(5), respectively. The line width value 
which inputted into the ETC is 9.4 km s -1 based fitting the shape of the 
line prohle using a Gaussian fitting method. We applied the Doppler shift 
correction in the ETC of V; sr = —22.54 km s -1 , this accounts for the local 
standard of rest radial velocity at the time of observation. Plots in Figure 8 
show that the flux calibrations of the emission lines from the ETC and the 
observational data are approximate. 

5. SUMMARY 

We have developed the IGRINS ETC to generate realistic simulations of 
both continuum and emission line sources. The signal and noise values are es¬ 
timated by taking into account telluric background emission and absorption, 
the emission and transmission of the telescope and instrument optics, and the 
dark current and read noise of the infrared detector arrays. We estimate the 
atmospheric transmission using models that account for the amount of pre- 
cipitable water vapor (PWV) along the line of sight to the target. From the 
ETC, the observers can estimate the S/N of the spectrum for each spectral 
resolution element given the exposure-time and number of exposures. 


The comparisons of the S/N vs. wavelength of HD 124683 and GSS 32 
show that the simulated S/N values from the ETC are overestimated by 
40—50%. Note that the adopted throughput values in Table 1 are mostly 
from rough estimates in ideal cases. In this paper, the comparisons of the 
simulated S/N values with that from the observed IGRINS data test the 
throughput dependent aspects of the ETC. 

The source-code of the IGRINS ETC version 3.10 is available to be down¬ 
loaded from the website http://irlab.khu.ac.kr/~igrins. 
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Table 1: Parameters in throughput estimation 


r 

Parameter 

Value 

Tat mo 

Atmosphere® (H band) 

0.77 

T"atmo 

Atmosphere® (K band) 

0.83 

T~mirror 

Telescope Mirrors 

0.55 

r gr 

Immersion Grating 

0.88 

T"vph 

VPH Grating (H band) 

0.80 

Typh 

VPH Grating (K band) 

0.75 

^optics 

Other IGRINS Optics 6 

0.42 

Tslit 

Slit-Loss® 

0.64 

Pje 

Q.E. of Detector 

0.80 


a See section 2.1. 

b Window, a pupil stop, collimators, and camera lenses. 
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Table 2: Parameters in sensitivity simulations 


Parameter 

H band 

K band 

Unit 

Center Wavelength 

1.63 

2.22 

/iin 

Zodiacal Light a at latitude of 0° 

99 

73 

photon s -1 m -2 arcsec -2 

Zodiacal Light 3, at latitude of 30° 

50 

37 

photon s -1 m -2 arcsec -2 

Zodiacal Light 3 at latitude of 60° 

42 

31 

photon s' 1 m -2 arcsec -2 

Zodiacal Light 3 at latitude of 90° 

30 

22 

photon s -1 m -2 arcsec -2 

Half Moon light b 

54 

30 

photon s' 1 m -2 arcsec -2 

Full Moon light b 

595 

277 

photon s -1 m -2 arcsec -2 

OH light c 

3390 

1096 

photon s -1 m -2 arcsec -2 

Ambient Temperature 

275 

275 

K 

Detector Read Noise 

5 

5 

electron 

Detector Dark Current 

0.02 

0.02 

electron s -1 


a See section 2.2. We assume the heliocentric ecliptic longitude is d=90°. 
b See section 2.2. We assume the separation between target and Moon is 60°, and the altitude of the 

Moon is 45°. 

c See section 2.2. The OH value use in equation (5) is total OH line intensity in H and K bands. 
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Table 3: Average intensity values in Flat image 


Detector 

Order regions 

Inter-order regions 

Ratio 

H band 

4208 

698 

0.20 a 

K band 

2345 

650 

0.38 a 


a We calculate the ratio for both H and K bands by [(698+650)/2]/[(4208+2345)-(698+650)/2] =0.1. 
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Table 4: Standard Stars 


Standard star 

Observing date (UT) 

^exp 

rriK a 

Type 0 

W/ a 



(s) 

(mag) 


(K) 

HD 124683 

2014 May 26 

120 

5.552 

AOV 

9500 K 

HD 155379 

2014 Jul 12 

120 

6.520 

AOV 

9500 K 


a http://simbad.u-strasbg.fr/simbad/sim-fbasic. 
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COLLIMATOR 
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Figure 1: IGRINS optical design sketch. 











































Figure 2: Average values of the Moon light at different Moon-target separa¬ 
tions in the case of the half Moon phase (solid and dashed lines) and in the 
case of the full Moon phase (dash-dot and dotted lines) in H and K bands, 
respectively. 
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X-Position (pixels) 

Figure 3: Flat image in the H band (top plot) and K band (bottom plot). 
The red lines show the example of an order region in the echellogram. The 
inter-order regions are the remaining parts of the image. 
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Limiting Magnitude 


Figure 4: Plot of S/N vs. limiting magnitude for continuum point sources 
in the H band and K band with exposure-time per frame, t = 600 s, and 
number of exposures, N = 6. The seeing is 1.2 arcsec and the PWV is 2 mm. 
The Moon light is not included in the calculation. 


19 




















































Figure 5: Plots of S/N vs. wavelength in the H band and K band. We show 
a point blackbody (T = 6000 K) source with 12.1 mag in the H band and 
12.0 mag in the K band. The seeing is 1.2 arcsec and the PWV is 2 mm. 
The exposure-time per frame is 600 s, and the number of exposures is 6. The 
Moon light is not included in the calculation. 
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Figure 6: Plots of S/N vs. wavelength of HD 124683 from the observational 
data (top plot) and simulation data from the ETC (bottom plot). The sim¬ 
ulated data from the ETC assumes exposure times of t = 120 s, number 
of exposures N = 4, the K-magnitude m k = 5.552 mag, and the effective 
temperature of the source T = 9500 K. The seeing is 0.9 arcsec. The PWV 
is assumed to be 2 mm. The Moon light is not included in the calculation. 
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Figure 7: Plots of S/N vs. wavelength of GSS 32 from the observational data 
(top plot) and simulation data from the ETC (bottom plot). The simulated 
data from the ETC assumes exposure times of t = 240 s, number of exposures, 
N = 4, the K-magnitude, m^ = 7.324 mag, and the effective temperature of 
the source, T = 1360 K. The seeing is 0.8 arcsec. The PWV is assumed to 
be 2 mm. The Moon light is not included in the calculation. 
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Figure 8: Plots of observational data and the ETC simulation of NGC 7023. 
The left plots show the H 2 1-0 S(l) emission line, and the right plots show 
the H 2 7-0 0(5) emission line. We applied the Doppler shift corrections in 
the ETC, V i sr = -22.54 km s -1 , the local standard of rest velocity. The flux 
calibration is done by comparing to the signal from HD 155379. The Moon 
light is not included in the calculation. 
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Appendix A. MODEL ATMOSPHERE 


We consider telluric absorption and emission lines in the IGRINS ETC 
using a model to calculate synthetic telluric spectra (Seifahrt, Kufl & Zngl et 
al., 2010). The model spectra cover from 1.4—2.5 //in with precipitable water 
vapor (PWV) of 2, 4, and 8 mm for the weather conditions at McDonald 
Observatory. The original transmission atmosphere spectra are convolved to 
be the same spectral resolution as IGRINS, R = 40,000 in the H and K 
bands. Figure A.9 and A.10 show the atmospheric transmission in the H 
and K bands. 

We compiled lists of 249 bright OH emission lines from the literatures 
(Rousselot, Lidman & Moreels et ah, 2000; Oliva & Origlia, 1992). The 
original OH emission line data are convolved to be the resolution of R — 
40,000. Figure A. 11 shows the spectra of OH emission lines in the H band 
and K band. 
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Figure A.9: Plots of atmospheric transmission in the H band. Precipitable water vapor (PWV) of 2 mm (red 
lines), 4 mm (blue lines), and 8 mm (green lines) for the weather conditions in the McDonald Observatory. 




















































































































Atmospheric transmission 





2.105 2.110 2.115 2.120 2.125 2.130 2.135 2.140 2.145 2.150 2.155 2.160 2.165 2.170 2.175 2.180 2.185 2.190 2.195 2.200 


1.0 

0.8 

0.6 

0.4 

0.2 

0.0 




——'yryy-— 




.. I .... I .... I .... I .... I ... I I .... I ... I I .... I .... I .... I .. I . I I ... I . I .. I .... I ... I I .... I ... I I I ... I I I I 

2.205 2.210 2.215 2.220 2.225 2.230 2.235 2.240 2.245 2.250 2.255 2.260 2.265 2.270 2.275 2.280 2.285 2.290 2.295 2.300 



0.0 ^ 


^00 2.305 2.310 2.315 2.320 2.325 2.330 2.335 2.340 2.345 2.350 2.355 2.360 2.365 2.370 2.375 2.380 2.385 2.390 2.395 2.400 



2.400 2.405 2.410 2.415 2.420 2.425 2.430 2.435 2.440 2.445 2.450 2.455 2.460 2.465 2.470 2.475 2.480 2.485 2.490 2.495 2.500 


Wavelength (urn) 


Figure A.10: Atmospheric transmission in the K band. Precipitable water vapor (PWV) of 2 mm (red 
lines), 4 mm (blue lines), and 8 mm (green lines) for the weather conditions in the McDonald Observatory. 















































































































Figure A.11: Plots of OH emission lines in the H band (top plot) and K 
band (bottom plot). 
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